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ABSTRACT 

The nearby isolated neutron stars are a group of seven relatively slowly rotating neutron stars that 
show thermal X-ray spectra, most with broad absorption features. They are interesting both because 
they may allow one to determine fundamental neutron-star properties by modeling their spectra, 
and because they appear to be a large fraction of the overall neutron-star population. Here, we 
describe a series of XMM-Newton observations of the nearby isolated neutron star RX J0806.4— 4123, 
taken as part of larger program of timing studies. From these, we limit the spin-down rate to v = 
(-4.3 ± 2.3) X 10-1*' Hz s-^ This constrains the dipole magnetic field to be < 3.7 x lO^^ G at 2cr, 
significantly less than the field of ~ 10^'' G implied by simple models for the X-ray absorption found at 
0.45 keV. We confirm that the spectrum is thermal and stable (to within a few percent), but find that 
the 0.45 keV absorption feature is broader and more complex than previously thought. Considering 
the population of isolated neutron stars, we find that magnetic field decay from an initial field of 
<3 X 10^* G accounts most naturally for their timing and spectral properties, both qualitatively and 
in the context of the models for field decay of Pons and collaborators. 

Subject headings: magnetic fields — stars: individual (RX J0806. 4-4123) — stars: neutron — X-rays: 
stars 



1. INTRODUCTION 

Th e so -called isolate d neutron stars (INS; see iHaberll 
120071 and iKaplanI 120081 for reviews) are a group of seven 
(confirmed) nearby (< 1 kpc) neutron stars with low 
(~ 10^^ ergs~^) X-ray luminosities and long (3-11 s) spin 
periods. They stand out from the normal neutron star 
population because of their timing properties (which 
should not have influenced their met hod of discovery 
via soft X-ray emission, although see IHevl fc KulkaTa 
Il998[ ). The X-ray luminosities are consistent with cool- 
ing neutron stars of ages 0.5 Myr (jPage et al.l l2004) . in 



the cooling rate a nd thermal content of the neutro n star 
atmosphere fe.g. JPons. Miralles. fc GeppertI 120091 . and 
references therein). 

Dipolar magnetic field strengths can be esti- 
mated from coherent timing solutions, and we 
used X-ray observations to derive such solu- 
tions for four IN S, finding magnetic fields of 
1 - 3) X 10^^ G (iKaplan fc van KerkwiikI l2005allbl; 



van Kerkwijk fc KaplanI 2008| _ Kaplan fc van KerkwiikI 



rough agreement with their kinemat ic ages ( Waltc^ 20011: 
Kaplan, van Kerkwijk, fc Anderson 2002g: iMotch et al.l 
2005; Kaplan ct al. 2007; Motch et al. 2003). 

A crucial unknown in understanding the INS is their 
magnetic field strengths. This is imp ortant in un- 



dersta nding the X-ray spectra (Zanc. T uroUa. fc Draki 
20041: IMotch. Zavhn. fc Haberl 2003' |Ho et al.l 1200 



van Ke rkwijk fc Kaplan 2007; Hab erl 120071 ) and their 
thermal history (Hcyl fc Kulkarni) 119981 ). The X-ray 
spectra of the INS appear thermal, with temperatures 
of ~ 10*5 K and with, in aU but one source, broad absorp- 
tion features at energies of 0.2 to 0.75 keV. These features 
may give clues to the chemical composition and struc- 
ture of the INS surfaces, but also depend strongly on the 
magnetic field. For a hydrogen atmosphere the simplest 
possibilities are electron and proton cyclotron resonances 
or transitions between bound states of neutral hydrogen, 
all of which depend on the field strength. Even with dif- 
ferent compositions (or states) the field strength is still 
relevant. At the same time, the magnetic field can affect 
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20091 hereafter IKvK05aL b IvKKOal IKvK09l: also see 
van Kerkwijk et al.l 1200^! ). Here, we constrain the 
spin-down rate and hence magnetic field strength of 
the INS RX J0806.4-4123 (hereafter RX J0806) with 
a series of dedicated XMM-Newton observations. We 
also present a preliminary spectral analysis, but defer a 
detailed phase-resolved analysis to a later paper. 

RX J0806 was identified as a possible neutron star by 
iHaberl. Motch. fc PietschI ([1998) on the basis of a soft 
thermal spectrum an d the absence of an op tical coun- 
terpart. Using XMM. IHaberl fc ZavhnI ()2002[ ) confirmed 
that the spectrum was soft and blackbody-like and iden- 
tified a candidate 11.37 s periodicity. Further observa- 
tions confirmed this periodicity and also suggested that 
the X-ray spectrum was not purely a black body, but had 
abroad absorpti on feat ure at «0.45keV ([Haberl et al.l 
I2004bl: hereafter [H±01). In what follows, we assume a 
distan ce of 250 pc to RX J0806, derived bv lPosselt et al.l 
l[20?)l from a comparison of the X-ray absorption col- 
umn density with a model for the interstellar medium 
(using A^H = 1.0 X 10^°cm~^, consistent with what we 
infer) . 

With RX J0806, we now have spin-down measurements 
or constraints for five of the seven confirmed INS. Exam- 
ining the global properties of the INS relative to other 
related neutron star populations, we attempt to under- 
stand how the INS fit with the other groups. 
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TABLE 1 

Log of Observations and Times of Arrival 



Rev. 



Date 



Exp.'' 
(ks) 



Counts^ 



TOA^ 
(MJD) 



168 
618 
1542 
1544 
1551 
1562 
1621 
1631 
1633 
1649 
1705 
1705 
1710 



2000 
2003 
2008 
2008 
2008 
2008 
2008 
2008 
2008 
2008 
2009 
2009 
2009 



Nov 08 
Apr 24 
May 11 
May 15 
May 29 
Jun 20 
Oct 15 
Nov 04 
Nov 09 
Dec 10 
Mar 31 
Apr 01 
Apr 11 



15.6 

22.0 
9.0 

10.0 
5.0 

13.0 
9.0 
5.0 
6.0 
9.0 
9.0 

26.0 
8.0 



25,604 
37,300 
11,467 
12,590 
7,107 
18,016 
12,224 
7,272 
9,386 
12,427 
11,185 
41,755 
11,027 



1.6 
1.5 
4.1 
4.4 
4.8 
9.2 
4.2 
4.7 
27.8 
4.4 
4.6 
23.7 
5.8 



51856 
52753 
54597 
54601 
54615 
54637 
54754 
54774 
54779 
54810 
54921 
54922 
54932 



,691599(3 
,747077(3' 
,501552(4 
,309174(4 
,274494(6 
,594919(5 
,488408(4' 
185892(8 
,285979(6 
,460578(4 
,908923(4 
,901462(4 
,057781(4 



Note. — All observations used the small window mode and thin 
filter for both EPIC-pn and EPIC-MOSl/2, except for Revs. 168 and 
618, in which the full window mode was used, which meant that only 
the EPIC-pn data were suitable for timing. 

^ The exposure time, number of counts, and estimated fraction of 
events due to background /bg given here are for EPIC-pn only. 
^ The TOA is defined as the time of maximum light of the funda- 
mental (following Eqn. [TJ closest to the middle of each observation 
computed from the combined EPIC-pn and EPIC-MOSl/2 datasets, 
and is given with l-cr uncertainties. 

The structure of this paper is as follows. In Section [2l 
we present our new data and perform our timing fij 12. ip 
and spectroscopic 12. 2|) analyses. We then consider the 
implications of these results. In § 13.11 we compare the 
timing and spectral properties of RX J0806 to the rest 
of the INS and try to understand the origin of the X-ray 
absorption features. In § I3.2[ we extend the compari- 
son to the pulsar population as a whole, exploring the 
qualitative relation between the INS and relevant sub- 
populations of pulsars. From this, we are led to consider 
the role that magnetic field decay may play in the evo- 
lution of the INS, and in ^ 13.31 we discuss our results in 
the context of a specific model for field decay, that of 
iPons et al.l (|2009f ). Finally, we conclude in §|4l 

2. OBSERVATIONS & ANALYSIS 

We observed RX J0806 eleven times with XMM 
pansen et all 120011) in 2008 and 2009, and focus here 
on the data taken with the European Photon Imaging 
Camera (EPIC) with pn and MOS detectors, all used 
in small window mode with thin filter s (Table HI). All 
our ob serv ations, as well as those from lHaberl &: ZavlinI 
(|2002f) and lH+041 (taken with the same filter, but with 
the full window mode instead) , were processed with S AS 
version 8.0.1. We used epchain and emchain and se- 
lected source events from a circular region of 37'.'5 radius 
with energies between 150 eV (the default) and 1.2 keV 
(where flares are minimized; the source is not detected 
above 1.2 keV, and the background grows increasingly 
dominant below 150 eV). As recommended, we included 
only one and two-pixel (single and double patterns 0-4) 
events with no warning flags for pn, and single, double, 
and triple events (patterns 0-12) with the default flag 
mask for MOSl/2. We barycentered the event ti mes us- 
ing the Chandra X-ray Observatory position from lH-l-041 : 
a = 08'^06"^23M0 and 5 = -41°22'30'.'9 (J2000). We 
extracted background lightcurves for pn from similarly 
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Fig. 1. — Pulse profile of RX J0806, repeated twice for clarity. We 
show the background-subtracted EPIC count-rate in each of 16 bins 
as a function of pulse phase. The solid curve is the best-fit model 
including the first harmonic, while for comparison we also draw the 
best-fit sinusoid as the dashed curve. The histogram/points include 
all events from our 2008—2009 observations phased up following 
Table [2] while the open circles are the longest single observation 
(Rev. 618 from 2003). 

sized regions offset from the source, but at the same RAWY 
coordinate, as recommended by the SAS User Guide. ^ 
For MOSl/2, the small-window mode does not permit 
such large background areas, but we used several smaller 
areas to compensate. The background rate for our data 
(up to ~ 24% in one case) was higher than for the archival 
data (^^2%), in some cases due to severe flares. However, 
restricting the time ranges to remove the flares did not 
significantly alter our results or precision. 

2.1. Timing Analysis 

As a starting place, we first determined the frequency 
that maximized the power in a periodogram for the 
EPIC-pn data from Rev. 1562 (the observation in 2008- 
2009 with the highest Zi power). We found a best- 
fit frequency oi v = 0.0879 58 ± 0.000006 Hz, consistent 
with that found by lH-l-041 for the earlier data. A Z| 
periodogram that incorporates the first harmonic (see 
below) gives the same result with slightly better pre- 
cision, where the uncertainty on v is found with the 
same analytic expression as for the Z"^ res ult (based on 
IRansom. Eikenberrv. fc Middleditehl I2002f) and we have 
verified with simulations that this is correct. 

Using the above frequency, we determined the times-of- 
arrival (TOAs; see Table [1]) for the combined EPIC data 
from each observat ion by fitting the binned lightcurves 
(following KvK05b: we verified that we obtained similar 
results from unb i nned fits following ICashI 119791 ) . Like 
IHaberl fc Za'^ (l200l) . we found that the lightcurve of 
RX J0806 was best described not by a single sinusoid but 
instead by a sinusoid and its first harmonic, where the 
number of counts in bin i are: 

N,^A {cos [27r(/oi, - 0o)] + ^2 cos [47r(/ot. - 0o - M2)]}+C 

(1) 

where /o is the frequency of the fundamental from above, 
ti is the time of bin i on the interval [0, l//o), and the 
fit parameters are amplitude A, phase (j)Q, amplitude ra- 
tio r2, ph ase offset and constant rate C, similar to 

lKvK05bl We found initial best-fit parameters from the 

See jhttp : //xmm.esac ■ esa. int/external/ximn_user_support/documentation/ 
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Fig. 2. — Phase residuals for RX J0806. We show the residuals 
relative to a linear model (i> = 0). The line shows the best-fit 
quadratic solution. 

longer observations, finding that all were reasonably con- 
sistent with a constant pulse profile with A/C « 4.6%, 
r2 « 0.45, and A(/)2 ~ 0.11 (corrected for background). 
T his is c onsistent with the ~ 6% pulsed fraction found 
bv lH-l-O-l . We refined the parameters using the combined 
event list derived from all 2008-2009 observations phased 
together (see below). 

The spacing and precision of the TOAs is insuffi- 
cient for an unambiguous timing solution (unlike in 
IKvKOSal b but hke in ivKKOS, and .KvK09ft . Instead we 
searched for possible coherent timing solutions by it- 
eratively trying sets of cycle counts between the dif- 
ferent TOAs fro m 2008 a nd 2009 (similar to ,vKKO S, 
although as with lKvK09l we did not incorporate fre- 
quency information from each observation because it 
does not add extra information). We limit solutions to 
|i>| 5, 1-5 X 10~^^Hzs~^, the 2>-a incoherent limit from 
IH+M . 

We find one solution that is considerably better than 
the alternatives with = 5.1 for 8 degrees of freedom 
(DOF) and a small, not-quite significant spin-down rate 
of = (-4.3 ± 2.3) X 10-16 Hz s'^ (see Fig. H Table O. 
This gives a 2 cr limit on the magnetic field of field of 
-Bdip < 3.7 X 10^3 G. The next best solution has = 27.2 
and a very different spin-down rate: (-1-4.32 ± 0.02) x 
10-i'*Hzs^i. For a fit with 3 free parameters, a change 
in of 22.1 means that the best-fit solution is favored at 
99.994% confidence. Even less likely solutions are found 
for combinations of other cycle counts. 

We confirmed this solution using a coherent Zf peri- 
odogram as a function of both u and i>. The best solu- 
tion was consistent with that found in the TOA analy- 
sis, although it varied by ^ 1 cr in 0. This is likely be- 
cause the Zl did not incorporate any information about 
the harmonic, while the TOA analysis did. The best-fit 
peak has = 162.3, consistent with the background- 
corrected rms pulsed fraction of « 3.8%. Spin-down is 
not detected with Zi{i> = 0) = 162.2, as the Zf does not 
include the harmonic. Including the harmonic at arbi- 
trary phase and amplitude we get Z| = 206.9 at {i/, v) = 
(0.0879477628±8x lO-i^'Hz, (-2.0±3.6) x lO'^^ Hz s'^), 
very close to the location of the coherent solution (the in- 
creased uncertainty comes from not including constraints 
on the pulse profile). 

After performing our coherent timing solution, we 



TABLE 2 
Measured and Derived Timing 
Parameters for RX J0806. 4-4123 



Quantity 


Value 




o^,oyo— o^,yoz 


to (MJD) 


54771.319605(2) 


V (Hz) 


0.0879477624(9) 


V (10-1^ Hzs-i) 


-4.3(23) 


TOA rms (s) ... 


0.4 




5.1/8 


P¥) 


11.37038592(12) 


P (10-"ss-i). 


5.5(30) 


T"char (Myr) 


3.3 


Sdip (10" G) . . 


2.5 


E (103°ergs-l) 


1.5 



Note. — Quantities in parenthe- 
ses are the formal l-tr uncertainties 
on the last digit. Tchar = P/'iP is 
the characteristic age, assuming an 
initial spin period Po "C P and a con- 
stant magnetic field; Bj;p = 3.2 X 

10^^ V PP G is the magnetic field in- 
ferred assuming spin-down by dipole 
radiation; E = 3.9 X 10*^ v!> erg s~^ 
is the spin-down luminosity. In addi- 
tion to the nominal values above, the 
2-a limits to those quantities (based 
on i> > -8.9 X lO-i'^Hzs-i) are 
Tchar > 1.5 Myr, Bdip < 3.7 X lO" G, 
and B < 2.8 X lO^Oergs"!. 

made a pulse profile using the events from all 2008- 
2009 observations (Fig.[T]). We found that the shape was 
consistent with what we had assumed above, but that 
changing the pulse shape parameters within the mea- 
sured uncertainties changed the timing solution slightly, 
typically by 10% in v (20% of the uncertainty on v). 
Similarly, changing the timing solution within the un- 
certainties changed the fitted parameters slightly. How- 
ever, since all such changes were significantly less than 
the uncertainties, we decided to choose an average set of 
pulse profile parameters and use them throughout. The 
final parameters that we used were A/C = 0.046 ±0.004, 
r2 = 0.45 ± 0.10, and A(/)2 = 0.11 ± 0.02 (i.e., our arrival 
times are based on fits in which only (pQ and C are free 
parameters). These values are consistent with the tim- 
ing solution given above and in Table [51 The profile with 
the harmonic is preferred greatly to just a sinusoid, with 
= 43.8 for 13 DOF (fundamental only) vs. = 14-2 
for 11 DOF (with the harmonic). Changing parameters 
slightly gives similar but slightly different results. For ex- 
ample, with r2 — 0.55 instead of 0.45, we find v decreases 
by 2% which is 5% of the uncertainty on v. Decreasing 
A/C to 4% increases v by 7% (15% of the uncertainty on 
v). The energy range that we used also affects our results 
slightly, but the solution remains consistent overall. For 
instance, we see slight differences if we restrict our anal- 
ysis to the energy range of 300-500 eV, where the pulsa- 
tions from RX J0806 are the very strong and the spec- 
trum shows X-ray absorption (see below) . While the pul- 
sation amplitude increases to A/C = 7%, the timing so- 
lution is consistent with i' = (-2.1 ±2.2) x lO^^^Hzs^^ 
Unfortunately, we cannot unambiguously extrapolate 
our solution back to the 2000 and 2003 observations and 
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Fig. 3. — Binned EPIC-pn spectra of RX J0806, where we have 
subtracted off the mean spectrum and divided by the uncertainties. 
The background has been subtracted, but no other corrections have 
been done. The individual observations are labeled by their dates 
from Table [1] 

thus infer a precise spin-down rate: the gaps are so large 
that even at 1 a, the uncertainty on the cycle count from 
the ly uncertainty is ^ai.At'^ = 7.2 cycles (where At = 
1.8 X 10® s is the gap between the reference time and 
Rev. 618). 



With 



2.2. Spectroscopic Analysis 

' 3 times longer total exposure time compared 
to H+04| . we wished to see whether the basic spectral fits 
of lH-l-04l are still valid and to look for possible long-term 
variability such as th at found for RX J0720.4— 3125 by 
Ide Vries et all (l2Q0l . 

To do so we examined all EPIC-pn spectra of 
RX J0806. (A full spectral analysis, including the EPIC- 
MOS and RGS data and a phase-resolved analysis, is 
in progress.) We used the same source and background 
extraction regions as for the timing analysis, created ap- 
propriate response files, and binned the spectral files such 
that the bin width was at least 25 eV (about one third 
of the spectral resolution) and the number of source plus 
background counts was at least 25. 

We first compared the raw EPIC-pn spectra of all of 
the observations against each other. This did not include 
any response files or calibration corrections, but even so 
the binned pn spectra were generally consistent with each 
other implying no spectral change (Fig. [3]). The only 
small deviations were seen for the archival data in 2000 
and 2003 (reduced of 2-3), and much of that could be 
corrected by a small slope across the 0.15-1.2 keV band: 
a fractional change of ~ 0.15 keV"""^ is sufficient to get 
a reduced of 1-2 for all observations. This change 
for the early data most likely refiects a small degree of 
event pileup present in those observations above 0.7 keV, 
caused by the 73 ms frame time of the full frame mode 
being insufficient to separate all events (for the small 
window mode, with its 6 ms frame time, pileup is not an 
issue). Because of this corruption, in the following fits 
we excluded the full frame data. We note, however, that 
the results including those data were still quite similar. 

We fit the pn da ta from 2008-2009 together using 
sherpa. Like lH-l-041 . we found that an absorbed black- 
body did not provide a good fit. Indeed, the deviations 
from a blackbody forced the interstellar absorption col- 
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Fig. 4.— EPIC-pn spectra of RX J0806. The top panel shows 
the data from each observation (as in Fig. [Sjl along with the best- 
fit blackbody plus 2 absorption line model (2G; solid line). The 
lower panels give the ratio of the data from 2008—2009 (with all 
observations averaged together for clarify as the blue points, and 
Rev. 618 separately as the red circles) to the model for the best- 
fit blackbody, IG, and 2G models (as labeled; Table [Sj. In the 
blackbody panel we plot the ratio of the models IG/Blackbody 
(dashed green line) and 2G/Blackbody (solid green line). Where 
we have Gaussian absorption components, we also show the Gaus- 
sians schematically (black lines), although the amplitudes are on 
an arbitrary scale. 

umn density iVu to zero. Rather than setting A'h arbi- 
trarily, we fit for an unabsorbed blackbody: while not 
realistic, this provides a reference. We list the corre- 
sponding radius and temperature in Table [3l In Fig. IH 
one sees that there are strong negative residuals near 
0.45 keV. We tried fitting for this using Gaussian absorp- 
tion lines, defined as follows, 



F{E) = Fc{E) 



1 - ^ Aicxp ( -41n2 



i=l 



' FWHMi^ 



(2) 

where Fc{E) is the continuum spectrum (a blackbody 
with interstellar absorption), Ei, Ai, and FWHMi are 
the central energy, absorption depth, and FWHM of ab- 
sorption component i and the factor of 4 In 2 converts 
from a FWHM to a Gaussian a. 

In Fig. m one sees that the absorption is broad and 
somewhat asymmetric. We found that, as a result, if we 
left the FWHM free, a single Gaussian became as wide as 
the blackbody peak. Instead, therefore, we decide d to fix 
the FWHM to 0.2 keV, similar to what was used bv lH-H04l 
(0.16keV). This is arbitrary, but should at least give a 
qualitative sense of the amount of absorption present. 
Including this single component results in a feature with 
an equivalent width of 53 eV and a significantly improved 
fit. As a consistency check, we kept the central energy 
and amplitude of the absorption fixed and fit for the 
FWHM, finding a best-fit value of 0.23 keV. 

In Figure 21 there still appear to be significant devi- 
ations of the residuals of the IG model near 0.65 keV, 
showing that our single line could not reproduce the de- 
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TABLE 3 

Results of Spectroscopic Fits to the Combined 2008-2009 EPIC-pn Data 



ModoP 




DOO b 










^2" 


Fx" 




X'/DOF 




(eV) 


(km) 


(1020 cm-2) 


(eV) 




(eV) 




(10-12 


ergs s ^ cm ^ ) 




Blackbody 


95.0(2) 


1.432(8) 


0.0 










2.3 


2.3 


649/392 


IG 


95.1(11) 


1.93(7) 


0.95(11) 


486(5) 


0.25(2) 






2.3 


2.9 


482/389 


2G 


87.2(11) 


2.39(15) 


1.7(2) 


460(5) 


0.42(3) 


693(12) 


0.26(3) 


2.3 


3.5 


402/387 



Note. — Only the 0.15-1.2 keV range was fitted. Quantities in parentiieses are tiie formal l-cr uncertainties on tlie last 
digit. Quantities without uncertainties were held fixed for that particular fit. 

^ A blackbody modified by interstellar absorption, plus 0, 1, or 2 Gaussian absorption lines following Eqn. [2] 

^ The effective temperature and radius of the best-fit blackbody as seen by a distant observer. The radius is scaled to a 

distance of 250 pc. 

The central energy and absorption depth of Gaussian absorption lines with fixed FWHM = 200 eV. The equivalent width 
is A(7r/41n2)i/2FWHM = 212.9AeV. 

Absorbed and unabsorbed fluxes in the 0.2-2.0 keV band. 



viations from the black body well. We therefore added 
a second Gaussian absorption component near that en- 
ergy.* This significantly improved the fit, with much 
flatter residuals and a reduced of 1-04. The change 
in is significant at the ~ lO"""^^ level, according to 
an F-test. We also fit Gaussians at a range of energies 
across the 0.2-1.1 keV band, and find that except for 
the edges of the band (where the fitting breaks down) 
only lines near 0.68 keV give significant improvements 
to the fit. Compared to allowing the FWHM of the 
first Gaussian to vary, we find that including a second 
line gives a much better fit (e.g., reduced of 1.52 for 
FWHMi = 0.5 keV). 

Even with two Gaussians, there is still structure in the 
residuals and the reduced is not quite 1, likely re- 
flecting the even-more complex shape of the absorption 
feature. The two lines in Figure [4] overlap, and there 
are residuals at their edges. The equivalent width of the 
first line has now increased to 89 eV, with the second at 
55 eV. The energies of these lines are in an almost 3:2 
ratio which may reflect some underlying harmonic struc- 
ture to the absorption, but it may also just come from at- 
tempts to fit the complex structure of the absorption with 
simple shapes: changing the FWHM of each line changes 
the central energies and the ratio, which is almost 1.7 for 
fines with FWHM of 0.3 keV. We could add additional 
components to improve the fit but they become increas- 
ingly artificial. We did try to fit in wavelength space 
wher e the shape of the absor ption is slightly different 
(e.g.. Ivan Kerkwiik et al.ll2004[ ). There the basic results 
hold, and while a single Gaussian does a somewhat bet- 
ter job, it is still far from perfect. Overall the absorption 
seems truly complex in shape, with a total equivalent 
width > 140 eV. Further phase-resolved fits may be able 
to help understand some of the complexities, but we defer 
them to a later paper. 

While the data were largely consistent with a constant 
spectrum, we did one simple test in which we allowed the 
individual observations to differ. With the basic 2G fit 
from Table [3] the observations had reduced ranging 
from 0.7 to 3.3, where the longest observation (Rev. 618) 
had the worst x^- This is also evident in Figured! where 

* IHaberll 1120071) discussed fitting two lines in a fixed ratio, and 
that improved the quality of the fit (albeit with less data) . How- 
ever, it was not well constrained, and both of lines from IHaberll 
(|20 07') would fit in the general region where wc find significant 
residuals. 



the residuals from Rev. 618 are generally similar to the 
average residuals, with some slight deviations at high 
energies. With the exceptions of the data from 2000 and 
2003, we find the data are consistent with a constant 
spectrum. The early data deviate slightly, giving slightly 
higher temperature and lower radius, but again this is 
largely because of the pileup-induced apparent hardening 
of the spectrum. Even including all of the data in the 
fit, kT changes by only a few eV with a resulting change 
in the blackbody radius such that -Rgg oc fcTj^g, i.e., the 
flux stayed relatively constant (< 2%). 

3. RX J0806. 4-4123 AND THE INS 

Much like the other INS for which we have derived 
timing solutions/constraints, RX J0806 has a magnetic 
field near a few by 10^^ G (nominally, 2.5 x lO^'^ G and 
< 3.7 X 10^^ G at 2 cr) with a characteristic age of a few 
Myr (nominally, 3.3 Myr and > 1.5 Myr at 2a) and a 
spin-down luminosity of ~ 10^° erg s~^ . As such, the INS 
form a very homogeneous group. Here, we discuss the 
implications for their spectral features, as well as for the 
overall neutron star population. 

3.1. INS spectral features 

In Ivan Kerkwijk fc KaplanI (|2007f l we compared the X- 
ray absorption of the different INS to a simple model in- 
cluding ion cyclotron and neutral hydrogen absorption 
in strong magnetic fields. We noted that the absorp- 
tion energy of RX J0806 («0.45kcV) was similar to that 
of RX J1605.3+3249 (van Kerkwijk ct al. 2004), and we 
see here that the depths are also comparable (we found 
139 eV for the equivalent wi dth of the mai n absorption 
line for RX J1605. 3-^3249). IHaberll (f2007l ) even found 
that a second line with an energy 1.5 times the first signif- 
icantly improved the fit of RX J1605.3-f 3249, as we found 
here for RX J0806, making the sources even more similar, 
although we again caution that it is not clear how unique 
or meaningful those fits are. We speculated that the mag- 
netic fields might be the same, at ~ 10"'^* G (this value is 
required for both p roton cyclotro n and neutral hydrogen 
transitions; also see lHaberlll200l . However, even though 
we do not have a clear detection of spin-down and thus 
cannot compare in detail, we can exclude this picture 
for RX J08()6: our limits show that at least the dipole 
component is quite a bit lower (see lKvK09l for a simi- 
lar discussion regarding RX J2143. 0-1-0654). The possi- 
ble presence of multiple lines complicates the situation 
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Fig. 5. — P-P diagram. We show the non-recycled pulsar popu- 
lation as points. The individual types of objects listed in Table |4] 
are also indicated: the INS are red circles (or upper limits where no 
P is known), the pulsars detected by ROSAT are blue squares, the 
nearby/young pulsars are blue five-pointed stars, and the pulsars 
with dipolc fields > 10^^ G and X-ray obse rvations are blue six- 
poin ted stars. We also plot the magnetars IWoods & ThompsonI 
1200 61) as green diamonds. We indicate lines of constant dipolar 
field and characteristic age as labele d. Finally , we also show the 
results of evolutionary model from Pons et al.l {2009) as the thick 
lines (solid lines: dipole only; dashed lines: including a toroidal 
component with a strength of 50 times that of the dipole). The 
initial dipole fields are I0l'''i3.5,i4,l4.5,l5 q ^j^g pQig^ there are 
cross-ticks at (0.01, 0.03, 0.1, 0.3, 1) Myr, we assumed initial peri- 
ods of 0.1s, and we have divided the field at the pole by 2 to get 
the field at t he equator for comparison with the spin-down field 
estimates fsee lLorimer &: Kramerll2004l ) . 

somewhat, since the "fundamental" may be at lower en- 
ergies and thus the inferred magnetic field weaker. How- 
ever, since harmonics should be significantly wea.ker fo r 
ion cyclotron resonance at least (jPavlov fc Panovlll976D . 
this seems unlikely to be the explanation. 

Presuming that the absorption comes from some tran- 
sition in the atmosphere, there are a number of alterna- 
tive explanations. For one, the magnetic field geometry 
could play a role. The dipole component that we mea- 
sure is only a projection of the true dipole field, and the 
field on the surface could have higher order o r substan- 
tial toroidal components (jBraithwaitd l2009f ). Beyond 
that, different INS could have different chemical com- 
position, although then we have to understand how the 
apparent emission r adii of the INS are relatively similar 
(jKaplan et al.ll2007D . 

Including RX J0806 in the magnetic field-effective 
temperature plane, it is consistent with the line de- 
fined by RX JI856. 5-3754 RX J0720.4-3I25, and 
RX JI308.6-I-2I27 (KvK09). As discussed in KvKOQ,, this 
correlation, while quite possibly a coincidence, may have 
some relation to the origin of the INS and the coupling 
of the magnetic and thermal evolutions (see below). It 
is also possible that it is influenced by surface condensa- 
tion^ (Mcdin & Lai 2007). 

3.2. The INS, and the Neutron Star Population 

While the INS were discovered over 10 years ago, we 
still lack a detailed appreciation for their place in the 

^ In our discussion of the effects of condensation in IKvK09l we 
incorrectly stated that a condensed surface would inhibit a vacuum 
gap and hence radio emission: in fact it is the opposite (Z. Medin 
2009, priv. comm.). 



overall neutron star p opulation and for wha t makes them 
unique (although see iPopov et all l2000bllal . [200l [20081 
and references therein for some discussion) . A full under- 
standing requires detailed models of the birth and evo- 
lution of neutron stars, but given the small numbers of 
objects and the large number of free parameters, getting 
reliable constraints is necessarily difficult. However, we 
can gain some insight by comparing the INS to different 
related neutron star populations. For this purpose, we 
give the salient properties of different objects in Table IH 
Aside from the seven confirmed INS, we consider: 

• Rotation-powered pulsars detected in the ROSAT 
All-sky Survey with count-rates > 0.05 s^^ in the 
PSPC (referred to here as "RASSPSRs"). 

• Other relatively young and nearby rotation- 
powered pulsars (referred to here as "NearPSRs" ; 
we adopted a distance limit of 1 kpc and a charac- 
teristic age limit of 5 Myr) . 

• Rotation-powered pulsars with dipole magnetic 
field > 10^3 Q (so called high-S pulsars, or "HBP- 
SRs"), but limited to those with X-ray observa- 
tions. 

These objects each make appropriate comparisons with 
the INS. The RASSPSRs are generally young and nearby 
neutron stars (we exclude the Crab pulsar from further 
consideration as it is much younger and more distant 
than the rest of the RASSPSRs), with X-ray luminosities 
and distances similar to those of the INS. Moreover, these 
were all detected in the same survey that discovered the 
INS. However, we do not know exactly how young the 
INS are. We therefore also include the NearPSRs: a 
sample of moderately young, moderately nearby pulsars 
regardless of their X-ray luminosity. While the INS are 
likely within 500 or 700 pc, we extend our range to 1 kpc 
for the pulsars since the distances of both classes are 
uncertain. Similarly, we extend to characteristic ages of 
5 Myr, compared to characteristic ages of ~ 3 Myr for 
the INS and kinematic ages of < 1 Myr. Finally, as the 
INS have larger than average magnetic fields, we also 
include those HBPSRs that have comparable dipole fields 
(> 10^'^ G) and X-ray observations. 
For each object, we give: 

• The timing properties, comprising the spin-period 
and the properties derived from that and P: 
dipole magnetic field, characteristic age, and 
spin-down luminosity. With the exception of 
RX J2143. 0+0654 and RX J0806 (where the spin- 
down measurements are only marginally signifi- 
cant) these data are uniformly of high quality. 

• The X-ray spectral properties, comprising the 
count-rate (as detected in the ROSAT PSPC), 
temperature and radius of the best-fit blackbody, 
luminosity of the blackbody, and total luminos- 
ity including any non-thermal components. The 
quality of these parameters varies by object. For 
the INS they are relatively uniform, but the de- 
tails of the spectral fit depend on the observa- 
tions and the object as different objects have dif- 
ferent levels of X-ray absorption (for example see 
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Ivan Kerkwiik fc KaplanI [20071 . §§ O andEU. For 
the RASSPSRs the thermal components are from 
just one possible decomposition, and often two 
blackbodies are required in addition to a non- 
thermal component. For the fainter NearPSRs and 
HBPSRs again there are problems of decomposi- 
tion (now limited by signal-to- noise ratio), and in 
many cases we only have upper limits. 

• The distance (measured by astrometry where 
available) and any other age indicator (kine- 
matic age or age of associated supernova rem- 
nant). Distances from astrometry are ideal, 

but in marginally-s ignificant ca ses (such as 

RX .T0720. 4-3125: iKaplan et al.l 120071) L utz- 
Kelker bias (|Lutz fc Kelkedll973HSmithll2003 fl can 
make them appear closer than they are, although 
other information such as X-ray absorption and 
Galactic geometry can be incorporated to improve 
the situation. For objects without astrometry, we 
rely on the dispersion measure (in radio) and hy- 
drogen column densities (in X-rays) along with 
models for the Galaxy. In both cases, the results 
can be unreliable, especially for close objects (as 
these are), although it is less likely that they are 
systematically biased. Furthermore, the X-ray ab- 
sorption column density can be unreliable as it is 
often covariant with other fitting parameters as 
well as with the assumed shape o f the spectrum 
(e.g., iDurant fc van Kerkwiilj|2006D . 

For additional notes on the data, see Table IH 

If we take the RASSPSRs and INS as a single sam- 
ple detected by ROSAT, the two groups have roughl y 
comparable sizes (IPopov et al.l[200a lKaplan"2004'. "2008^ . 
Based on their temperatures and independent ages at 
least some of the RASSPSRs may be slightly younger 
than the INS, but the temperatures can be affected by 
non-thermal emission processes absent in the INS. Com- 
paring the INS to the RASSPSRs, the INS h ave compara- 
ble X-ray lum inosities of 10"^^ ergs" ^ (e.g..lP opov et all 
[200llKaplanii2004l 2008). The kinematic ages of the INS 
(when available) are comparable to if not slightly larger 
than the ages of the RASSPSRs. However, the character- 
istic ages of the RASSPSRs are also small, < 1 Myr, while 
the characteristic ages of the INS are all > 2 Myr. All of 
the RASSPSRs have periods <0.4s, while the INS have 
periods an order of magnitude longer and considerably 
stronger magnetic fields.^ In contrast, for the slightly 
older NearPSRs the characteristic ages are often com- 
parable to those of the INS but the X-ray luminosities 
are several orders of magnitude lower lO^^ergs"^). 
The INS also have much lower E than any other popu- 
lation. Finally, we note that INS have larger radii and 
smaller effective temperatures for the same X-ray lumi- 
nosity. Is this an innate difference, perhaps as a result of 
the magnetic field geometry? Or is it artificial, perhaps 
a consequence of unrealistic emission models? 

® At some level, our comparison of magnetic fields, pulse peri- 
ods, characteristic ages, and spin-down luminosities is degenerate 
since these parameters are all related. Nonetheless, we will con- 
tinue to emphasize difli'erences among multiple parameters when 
appropriate. 



As we discussed previously (e.g., iKvKOSal) . we need 
to reconcile the discrepant kinematic and characteristic 
ages of the INS. For one or two objects it could be a 
coincidence, but even for the objects without kinematic 
ages (RX J0806, RX J2143.0-I-0654), the luminosities 
suggest the true ages are substantially shorter than the 
characteristic ages. While having a very long initial 
period is possible, it would have to be very close to the 
current value. For other sources where the character- 
istic ages exceed the true ages , like PSR J0538+2817 
(|Kramer et al.l l2003at iNg et all [20071 ). the initial spin 
period required in order to have spun-down to their 
current periods in their true ages (assuming dipole 
braking) is < 0.2 . T his is longer than the traditionally- 
assume d ^ 30 ms (iLyne. Pritchard. fc Graham-SmithI 
119931 iMighazzo et al.i l2002f ). but similar to the 
more recent ly favored pe riods o f a few hun- 
dred ms ([Kramer et al.l l2 003at iVranesevic et ahl 
200l ^ [G otthclf, Halpcri rfcSewardI [200a 
Faucher-Giguere fc Kaspi ,2006i) . and, most relevant 
for the present purposes, much smaller than what would 
be require d for the IN S (^5s). 

Kramer et al.l (|2006l ) have discovered a pulsar whose 
spin-down torque varies quasi-periodically along with its 
radio emission. For approximately 10-20% of the time 
PSR B193H-24 (an otherwise ordinary pulsar with char- 
acteristic age of 1.6 Myr and dipole field of 2.6 x 10^^ G) 
is visible as a radio pulsar, but for the remaining time 
the radio emission is not detectable and the torque is 
reduced to | of its normal value. The difference is at- 
tributed to the presence (or strength) of an energetic 
plasma wind, where a strong wind leads to radio emis- 
sion and higher torque. While we suspect from its Ha 
nebula that RX J 1856. 5— 3754 does indeed have an ener- 
getic wind and brakes by magnetic dipole radiation, that 
scenario is not completely self-consistent (see below). So 
the absence of radio emission from the INS may point to 
cessation (or at least diminution) of a wind and its asso- 
ciated torque, and this could be a recent development in 
the histories of the sources if radio emission has just shut 
off. If the wind/radio emission only stopped recently (the 
worst case scenario) then the period evolution would have 
been dominated by a spin-down approximately | times 
what we see now and the characteristic age that we mea- 
sure now would be high by 50%. This would go some 
but not all of the way toward resolving the discrepancy 
between the timing and ki nematic ages, although differ- 
ences in geometry (also see lKvKOgI ) could lead to slightly 
larger factors (other intermittent pulsars have since been 
discovered with slightly larger ratios of 1.7 instead of 
1.5, possibly due to differen ces in align ment between the 
rotation and magnetic axes; lKramej[2008t) . Such a com- 
parison also suggests that there could be torque and/or 
radio flux variations in the INS, neither of which has 
been seen (except for the one timing/spectra l change 
for RX J0720.4-3125: lvan Kerkwiik et al.|[200l but the 
sampling has been very sparse. 

Another possibility is that spin-down does not fol- 
low the expectations for magnetic dipole radiation. For 
instance, a decaying magnetic field would give rise to 
a situation like what we see for the INS, as it would 
allow the sources to spin down rapidly early in their 
evolution (|KvK09f) . Field decay was proposed previ- 
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ously for the INS as a way to keep them hotter longer 
and thereby make them ov errepresented in a local sam- 
ple ()Hevl fc Kulkarnilll998D . Initial timing results sug- 
gested that field d ecay was not p resently heating the INS 
(jZane et alJ 120021 : iKaplan et an i2002bl. and indeed the 
INS have roughly the expected thermal luminosities for 
their kinematic ages (especially given the large uncer- 
tainties and the steep decli ne in thermal lum inosity for 
objects of this age; see e.g. JPage et al.ll2^004[) . However, 
as we discuss below, field decay may have profoundly af- 
fected the magnetic and rotational evolution of the INS, 
and may have some other visible consequences. The 
models of field decay used bv IKaplan et al.l (|2002b[ ) were 
rudimentary, and including decay by other modes as well 
as a strong toroidal component can greatly change the 
outcome. 

Regardless of the specific model, if the characteris- 
tic ages are correct for the INS, they are too luminous 
by more than a factor of 100 to be powered by resid- 
ual heat. They would therefore require an extra energy 
source which we know cannot be the spin-down luminos- 
ity: the remaining alternative is magnetic field decay. If 
instead (and as we believe) the kinematic ages are cor- 
rect, then the X-ray luminosities are reasonable, but we 
need to explain the long periods (and the associated large 
characteristic ages, low spin-down luminosities). Again, 
magnetic field decay seems to be the best option. 

3.3. The INS and a Model of Magnetic Field Decay 

The coupled magnetic and therma l evolution f or neu - 
tron stars was studied in detail by iPons et all (|2009D . 
These authors assumed a range of initial field and tem- 
perature configurations and then followed them over 
time. Briefly, they found that for all sources with ini- 
tial magnetic fields >5 x 10^^ G, the final magnetic field 
is ~ 3 X 10^^ G at 0.5 Myr. But for all sources with weaker 
initial magnetic field the final magnetic field is j ust half 
the in itial field. We note that the model of IPons et al.l 
(|2009f) assumes magnetic field configurations and mech- 
anisms of decay that may well be overly simplified. One 
of the biggest free parameters in their models is the ratio 
of tor oidal to polo idal fields. This can vary significantly 
(e.g.. lBraithwaitel[200al . and the non-linearity of field 
decay means that a strong, decaying toroidal component 
can alter the weaker poloidal component while being oth- 
erwise invisible. 

Despite the above uncert ainties, i t seems encouraging 
that the main findings of Pons et al.l ([2009) for the results 
of field decay are similar to what we infer for the INS: 
neutron stars grouped around 2 x 10^^ G at true ages of 
~ 0.5 Myr or older. Taking the magnetic field evolution 
and using it to infer the spin-down history, one expects 
that neutron stars with strong fields quickly move across 
the P — P diagram, losing memory of the initial spin pe- 
riod. We find that as long as the initial field is between 
1 X 10^^ G and 7x 10^'' G, then the neutron star ends up at 
a period of 3 to 15 s (Fig. [5]). Coupling the thermal evo- 
lution in with the magnetic field. Pons et al. (2009) find 
luminosities of ^ 10'^^ erg s~^ at ages of ~ 0.5 Myr, similar 
to what is expected for less strongly magnetized neutron 
stars; the excess energy from the field decay is mostly 
radiated at earlier ages, with the thermal luminosity be- 
ing substantially higher around 10^ yr. At ^ 0.5 Myr, 



the difference may still be a factor of a few (dependent 
on field configuration and strength), but this would be 
largely lost in the observational uncertainties (age, dis- 
tance, etc.). This small remaining difference should not 
lead to a great over-representation in a local sample, and 
therefore the fact that about half of the young neutron 
stars detected by ROSAT (the INS plus the RASSPSRs) 
have long periods (and presumably strong fields) should 
give a reasonable clue to the true population. In a bit 
more detail, this will depend on the slope of the cool- 
ing curve: for power- law cooling L oc increasing the 
luminosity by a factor of A will lead to an increase in 
population size of A^^^" in a flux-limited sample if the 
population is in a constant volume (as might be expected 
for soft X-ray sources, since the exponential cutoff of in- 
terstellar absorption limits their detectability to < 1 kpc) . 
If the volume can increase too, we get an additional fac- 
tor of A in the population (assuming the population is 
confined to the Galactic disk). The true situation will 
probably be in between these extremes. With a typically 
between 2 and 3 at a few 10^ yr (for photon cooling), a 
factor of at most a few increase in luminosity (due to 
field decay) would lead to a factor of < 2 increase in pop- 
ulation. Therefore, the INS would represent between one 
quarter and half the total population. 

As an aside, we n ote that an interesting implica- 
tion of the model of IPons et all ()2009D is that neu- 
tron stars may be hotter on the equator than at the 
poles, a combination of Joule heating from field decay 
and reduced conductivity preventing the heat from go- 
ing inward. While the effect may be quite small at 
the ages of the INS, this may still complicate the in- 
terpretation of the lightcurve and phase-resolved spec- 
troscopy of the INS. Typically, a mo del with a warm 
pole and cool e quator is assumed fe.g.. lBraie fc Romanil 
l200llHoll200l . although more complicated models have 
also been considered (|Zane &: Turollal [20061 ). If the ge- 
ometry were reversed that could lead to different in- 
terpretations, and work is on-going to see if there are 
clear consequences to that (W. C. G. Ho, 2009, pers. 
comm.). The asymmetry between equator and pole 
is time- and field-dependent, being most apparent for 
strong fields and at early times. For weak fields, the 
standard hot poles are regained. It may be that the 
very low pulsed fractions and limits found in sources like 
RX J1856. 5-3754 and RX J1605.3-I-3249 are a conse- 
quence of being near the time where the surface becomes 
nearly isothermal, and that this underlies the difficultly 
in finding geometries that satisfy the observed pulsa- 
tion limits (|Braie fc Romanil 120021: IHoI|2007D . If there 
is still a strong, buried toroidal field it could also lead 
to some ongoing, low-level decay, which might give rise 
to some of the spe ctral and temporal evolution seen in 
RX J0720.4-3125 (lde\^ies et al.ll2004l: IVink^et al.1 120041: 
IHaberl et all 120061: Ivan Kerkwiik et al.ll2007f) . although 
nothing similar has been s een in the long-term monitor- 
ing of other sources f^ lHH IHabe rl 2007; KvK09). 

If, as we posit, magnetic field decay has influenced 
the period evolution of the INS, we can ask what 
the progenitors and de scendants of the INS might be. 
IHevl fc KulkaTml (fT998h discussed RX J0720.4-3125 as 
an old magnetar. Thi s could still be th e case, although 
based on the models of I Pons et al.l (|2009( ) the initial mag- 
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netic fields for the INS would be close to ~ 3 x 10^"* G 
and not much higher and would have followed an evo- 
lution that only skimmed the parameter range occupied 
by the magnetars (see Fig. [5]). The descendants of the 
stronger-field magnetars would be expected to end up 
with longer periods. The INS do not actually seem to be 
old versions of at least the X-ray-bright HBPSRs. These 
HBPSRs mostly have E > lO^'^ergs"^ (although this 
may partly be a selection effect, as this wo uld increase 
Lx)- Using the model of iPons et al.l (|2009[ ) and evolv- 
ing typical HBPSRs to ages of 0.5 Myr, we would expect 
E ^ lO^^ergs"^, which exceeds by an order of magni- 
tude what the INS have. These HBPSRs may in fact 
have magnetic fields close to their initial values with lit- 
tle decay, while the INS started with > 10^^ G and have 
decayed. A few of the HBPSRs, hkc PSR B0154+61 
[logioE = 32.8, logioLx < 32.3), PSR J1819-1458 
(logio E = 32.5, logio Lx = 33.7), and PSR J1718-3718 

(logj^o E — 33.3, logj^Q Lx ~ 33.5), may be more similar, 
although we have no independent age estimates for those 
objects. Of these, the comparison to PSR J1819— 1458 
may be particularly interesting, as it emits only sporadic 
radio bu rsts as a so-called "Rotating RAdio Transient" 
(RRAT; iMcLaughlinlFaLl [2006,1 . The possibihties of a 
connection between the RRATs and the INS have al- 
ready been discussed in several places fe.g.. 'Popo v et al. 
2006), but the spectral similarities (McLaughli n et al. 
20071 ) and the recent detection of extended X-ray emis- 

sion from PSR J 1819— 1458 that seems too large for its E 
(|Rea et al.ll2009D highlight it even further. The connec- 
tion can only go so far, though: searches for RRAT-like 
radio emission from th e INS have not been successful 
(jKondratiev et al.| [200?). 

As for the descendants of the INS, while they might 
be too faint for X-ray detection, they could still be ap- 
parent in radio surveys (although there are observational 
selection effects against long-period objects). However, 

very few objects are known in that part of the P — P 
diagram, and there are not enough to span the expected 
ages of up to 10 ''yr. Does this mean that the old INS 
have very narrow radio beams that miss the Earth, there- 
fore reducing their prevalence in radio surveys? Or are 
there no radio beams? Will the surfaces condense to al- 
low the formation of vacuum gaps (Medin & Lai 2007) 
or is this inhibited? This is largely a function of compo- 
sition and state, which, as we discussed, remains uncer- 
tain. The Ha nebula around RX J1856.5— 3754 might 
be taken as evidence for the generation of energetic 
particles (jvan Kerkwiik fc^ulk arni 2001, IKaplan et al.l 

I2002c[ ). but the requ ired E is fa r highe r than wh at is 
inferred from timing (jvan Kerkwiik fc Kaplanll2008f ). 

One interesting object that may bridge the gap 
between the standard pulsars and the INS is 
PSR B1845— 19: its long period and characteristic age 
put it withing reach of the INS, although we do not have 
an independent age or X-ray luminosity. Compared to 
the INS, the E for this object is slightly higher, while the 
characteristic age is similar, possibly suggesting that it 
sits at the boundary where magn etic field decay be comes 
important (> 10^^ G, according to lPons et al.ll2009l ). Just 
how close the actual age is to the characteristic age may 
prove an important test of this scenario. 



4. CONCLUSIONS 

We have presented initial results from coherent tim- 
ing of the nearby neutron star RX J0806. While we 
were not able to obtain a statistically significant mea- 
surement of spin-down, we were able to constrain the 
spin-down rate (and hence the dipolar magnetic field) to 
rather low values, and with the addition of a few data 
points with longer time baselines we should achieve a 
reliable measurement. The limit on the magnetic field 
(< 3.7 X lO^'^ G at 2 cr) is interesting, as it is considerably 
lower than the ~ lO^'' G expected from simple models of 
the X-ray spectrum . This e choes the discrepancy seen in 
RX J2143.0-t-0654 (lKvK09D . suggesting that we need to 
develop an improved model for the X-ray spectra of the 
INS. Whether this is just from an improved treatment of 
the radiative transitions (e.g., auto-ionizing transitions 
may play a role) or elements beyond hydrogen we do 
not know, but ongoing work in obtaining better timing 
constraints as well as phase-resolved spectroscopy should 
help to narrow the possibilities. 

The X-ray spectrum of RX J0806 shows clear signs 
of a broad absorption feature at 0.3-0.6 keV, and our 
data show that it cannot be simply modeled by one or 
even two Gaussians. The pulsations are also relatively 
strong in the same spectral region, giving some hope that 
phase-resolved spectroscopy will allow us to illuminate 
the possible decompositions. 

Comparing the INS with other relevant sub- 
populations of neutron stars, we are led to the conclusion 
that magnetic field decay has operated over the '^0.5 Myr 
lifetimes of the INS. This can be seen from a qualitative 
comparison of the observed X-ray luminosities of the INS 
with X-ray luminosities of other sources and with predic- 
tions from standard cooling curves. From both, one in- 
fers that the characteristic ages are systematically long, 
while the kinematic ages are consistent with cooling mod- 
els and the ages of other neutron stars with comparable 
luminosities. Th is conclusion is s upported by the de- 
tailed modeling of lPons et al.l ()2009D . whose expectations 
for neutron stars born with magnetic fields of ~ 2 xlO^* G 
evolved to ~ 0.5 Myr greatly resemble the INS. Among 
the implications of this model are that some neutron 
stars with magnetic field decay should be hotter at the 
equator than at the poles, in contrast with most assumed 
models of neutron star surface temperature distributions. 
Once again, phase-resolved spectroscopy would seem to 
be one of the best ways to try to discriminate between 
the possibilities. Improved models for the surface tem- 
perature, which constrain the viewing geometry, can then 
also be combined with proper motion measurements to 
try to understand some aspects of the complicated rela- 
tion between magnetic fields, rotation, and kick velocities 
in young neutron stars (Lai 2001). 

We thank A. Spitkovsky for helpful discussions. Based 
on observations obtained with XMM-Newton, an ESA 
science mission with instruments and contributions di- 
rectly funded by ESA Member States and NASA. DLK 
was supported by NASA through Hubble Fellowship 
grant #01207. 01-A awarded by the Space Telescope Sci- 
ence Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., for NASA, 
under contract NAS 5-26555. This research was sup- 
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TABLE 4 

Properties of the Isolated Neutron Stars and Related Rotation-Powered Pulsars 



Source 




Timing^ 








Sped 


rum^ 






A ffp'^ 


Refs. 




p 


-Bdip 


''"char 


-j 


PSPC 


kT 


^BB 


§10 ^X,BB 


iOgiQ ^X.tot 












(10" G) 


(Myr) 


(ergs-l) 


(s-i) 








V^^6^ ) 








Isolated Neutron Stars 


























RX J1308.6+2127 


10.31 


34 


1.5 


30.6 


0.3 


102 


4.1 


32.4 


32.4 


500 


0.8-1.2 


li 2i ^ ^ ^ ^ 


RX J0720.4-3125 


8.39 


24 


1.9 


30.7 


1.6 


87 


6.4 


32.5 


32.5 


360 


0.5-1.0 


7:'8r"9r"l0: ri~2 


RX J0806.4-4123° 


11.37 


25 


3.3 


30.2 


0.4 


92 


1.3 


31.2 


31.2 


250 




13; 14, this work 


RX J2143.0+0654 


9.44 


20 


3.7 


30.3 


0.2 


102 


3.2 


32.1 


32.1 


430 




15; 16; 17; 18; 19 


RX J1856.5-3754 


7.06 


15 


3.8 


30.5 


3.6 


62 


6.2 


31.9 


31.9 


160 


0.4 


20; 21; 22 


RX J1605.3+3249 










0.9 


93 


4.7 


32.3 


32.3 


390 


0.1-1.0 


23; 24; 14; 25; 26 


RX J0420.0-5022 


3.45 








0.1 


45 


3.3 


30.8 


30.8 


345 




i3;i4 


Rotation Powered Pulsars with PSPC > 0.05 s 


(RASSPSRs) 
















Crab 


0.03 


3.8 


0.0012 


38.7 


48.4 








37.6 


2000 


0.001 


27 


Vela 


0.09 


3.4 


0.011 


36.8 


3.4 


128 


2.1 


32.3 


32.8 


287 


0.01 


28; 29 


PSR B0656+14 


0.38 


4.7 


0.11 


34.6 


1.92 


56 


20.9 


32.7 


32.7 


288 


0.1 


30; 31; 32 


PSR B1951+32 


0.04 


4.9 


0.11 


36.6 


0.07 


130 


2.1 


32.2 


33.3 


2000 


0.06 


33; 34 


Geminga 


0.24 


1.6 


0.3 


34.5 


0.54 


43 


8.6 


31.9 


31.8 


250 




32; 35 


PSR B1055-52 


0.20 


1.1 


0.5 


34.5 


0.35 


68 


12.3 


32.6 


32.6 


750 




36; 32 


PSR J0538+2817 


0.14 


0.7 


0.6 


34.7 


0.06 


160 


3.2 


32.9 


32.9 


1470 


0.04 




Rotation Powered Pulsars with d < 1 kpc, r^hai 


< 5 Myr 


(NearPSRs) 














PSR J1741-2054 


0.41 


2.7 


0.4 


34.0 










31.1 


400 




39 


PSR B0450-18 


0.55 


1.8 


1.5 


33.1 












760 




40 


PSR B0450+55 


0.34 


0.9 


2.3 


33.4 












1190 




40 


PSR J1918+1541 


0.37 


1.0 


2.3 


33.3 












680 




PSR B0628-28 


1.24 


3.0 


2.8 


32.2 


0.003 






< 29.4 


30.1 


332 




41; 42 


PSR B2045-16 


1.96 


4.7 


2.8 


31.8 












950 


2 


40 


PSR B1845-19 


4.31 


10.1 


2.9 


31.1 












950 




PSR B0834+06'' 


1.27 


3.0 


3.0 


32.1 




170 


0.03 


29.0 


29.0 


640 




43 


PSR B1929+10 


0.23 


0.5 


3.1 


33.6 


0.012 


300 


0.03 


30.0 


30.4 


361 


1-2 


44; 45 


PSR J1908+0734 


0.21 


0.4 


4.1 


33.5 


< 0.005 










790 




46 


PSR B0823+26 


0.53 


1.0 


4.9 


32.7 


0.0016 


< 101 


Q 2g 


< 29.7 


29 4 


360 




^7 


PSR B0943+10'^ 


1.10 


2.0 


5.0 


32.0 




270 


0.02 


28.7 


28.7 


630 




48 


Pulsars With B > 10^^ 


G and X-ray Observations (HBPSRs) 
















PSR J1119-6127 


0.41 


41 


0.0016 


36.3 




210 


3.2 


33.4 


33.6 


8400 


< 0.01 


49; 50; 51; 52 


PSR B1509-58 


0.15 


15 


0.0016 


37.3 










34.3 


5200 


0.006-0.02 


53; 54; 5^ 


PSR J1846-0258' 


0.33 


49 


0.0007 


36.9 










34.5 


6000 


< 0.001 


56; CT,;^ 


PSR J1124-5916 


0.14 


10 


0.002 


37.0 




< 102 


12j 


< 33.3 


33.1 


6000 


0.03 


59; 60 


PSR J1930+1852 


0.14 


10 


0.002 


37.0 










34.6 


6200 




61; 62; 63 


PSR J1718-3718 


3.38 


74 


0.03 


33.3 




145 


7.7 


33.5 


33.5 


4500 




64 


PSR J1819-1458 


4.26 


50 


0.12 


32.5 




140 


11 


33.7 


33.7 


3600 




65 


PSR B0154+61 


2.35 


21 


0.2 


32.8 




< 63 


lOj 


< 32.3 


< 32.3 


1700 




66 



References. — [J: [Kaplan, Kulkarni, & van Kcrkwiik (2002a); 2: Kaplan & van Kcrkwiik (2005b): 3: Haberl ct al. (2003): 
[1 ISchwope et al.l |[2005l'); '5: SchwoDe et al. (2007); 6: Motch et al. (2009); 7: Kaplan et al. (2003b); 8: Motch ct al. (2003); 9: 
Kaplan k, ra n Kerkwijk (2005a); 10: Kaplan ct al. (2007); 11: van Kerkwiik ct al. (2007); 12: Haberl et al. (2004a); 13: Habcrl et alj 
Il2004bl); lit Habcrl (2007); 15: Zampieri et al. (2001); 16: Zane ct al. (2005); 17: Rea ct al. (2007); 18: Zane ct al. (2008); 1^ 
Kaplan &^an Kerkwiik (2009); 20: Burwitz et al. (2003); 21: van Kerkwiik & Kaplan (2007); 22: van Kerkwiik & Kaplan (2008); 
23; Kaplan, Kulkarni, & van Kerkwiik (2003a); 24: van Kerkwiik et al. (2004); 25: Motch ct al. (2005); 26: Zane et al. (2006); 27: 
Willingale et al. (2001); 28: Pavlov ct al. (2001); 29: Dodson et al. (2003); 30: Marshall & Schulz (2002); 31: Briskcn et al. (2003); 3^ 
De Luca ct al. (2005); 33: MigUazzo et al. (2002); 34: Li, Lu, & Li (2005); 35: Fahertv, Waher, & Anderson (2007); 36: Kramer et aQ 
T2p03b); 37: Roniani k, Ng (2003); 38: Nk ct al. (2007); 40: Chatterice et al. (2009); 43: Gil ct al. (2008); 44: Chattcriec ct al. (2004); 
g|: Misanovic. Pavlov, & Garmire (2008); 46: Becker & Triirnper (1997); 47: Becker et al. (2004); 48: Zhang, Sanwal, & Pavlov (2005); 
m-. Crawford ct al. (2001); 50: Caswell et al. (2004); 51: Gonzalez et al. (2005); 52: Safi-Harb & Kumar (2008); 53: Grcivcldinger et aO 
1^95); 54: Gacnslcr et al. (1999); 55: Gacnslcr et al. (2002); 56: Leahv & Tian (2008); 57: Gavriil ct al. (2008); 58 Kumar & Safi-Harfl 
j2&0a) ;|59i: iGaensler fc Wallacei (,2003i): j6Q: i Hughes et al., (i2003|): 1611: |Lu et al., (,200a): ,62i: ,Camilo et al.. (,200a ): ,63.: ,Leahv. Tian. fc Wani 
12008) ; IM: iKaspi k, McLauehlinI ll200g) T[65r TMcLauehlin et al.l ( 120071 ); |66|: IConzalez et al.l (120041) ; 

Note. — General pulsar data were taken from IPossenti et"an (120021 ) . iKaplan et al.l ( 120041) . and [Manchester et al.l (120051) . Within each 
class the objects are ordered by increasing characteristic age. 

Spin period P, dipole magnetic field -Bdip = 3.2 X 10^® \/ PP, characteristic age T^har = ^'/2P, and spin-down energy loss rate E = 
3.9 X 10""^ P/p3. 

^ ROSAT PSPC count rate (if known), effective temperature and radius (for spherical emission) of the best-fit blackbody component as 
measured at infinity, bolometric blackbody luminosity, and total luminosity in the 0.1-2.4 keV band. For sources with multiple blackbody 
components, we took the one with larger radius/smaller temperature. 

Parallax distance if known. If not, for the INS we use distances based on X-ray absorption (IPosselt et al.l 120071 ) or 500 pc for 
RX J1308. 64-2127, while for pulsars we use distances from dispersion measures dCordes fc Lazioll200g ). 

Independent age estimate. For the INS, this is the kinematic age, derived from tracing the object back to a probable birth site. The 
range is based on the observed range in distance as well as the result of multiple possible lairth sites. For pulsars this is either a kinematic 
age or an age estimate for the associated supernova remnant. 

Here we use the nominal values for the spin-down from Table [2] using the upper limits does not change the conclusions. 
^ The X-ray data are consistent with no non-thermal emission, but such emission cannot be ruled out. 
s The radius was assumed to be that of a polar cap. 

^ The emission from this source is also consistent with non-thermal emission with luminosity ~ 2 X 10^^ ergs^^. 

' The X-ray source is variable, and was observed at a level of about 10 X the previous flux (|Gavriil et al.|[2b08l ; IKumar k Safi-Har9l2008l ) . 
■> The radius was fixed to that of the whole surface. 
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